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ABSTRACT 


AM underwater acoustics™experiment conducted in shallow 
water (70 feet) off the New Zealand east coast in 1972-1973 
is described. Short acoustic pulses of 35 and 65 kHz sound 
Were projected along near-orthogonal paths of approximately 
300 yards. Environmental parameters were simultaneously 
observed. 

Statistical and spectral analyses of pulse heights were 
performed on 12 selected runs using digital techniques. 
Coefficients of Gaeta oiler ranged from 2.0% to 15.5%. Mai 
almost all cases, higher variability was observed along the 
acoustic path oriented perpendicular to the predominant 
Swell direction. Along this same path, periods corresponding 
to common surface swell periods were frequently evident in 
bhemaulocOrue late tons Tunctions of the filncluations= Coherence 
between the fluctuations along each path was low, averaging 
about 0.1. Long period oscillations suggestive of modulation 
by internal waves were apparent in several runs. No signif- 
icant dependence of variability on acoustic frequency was 
detected. 

Microscale temperature fluctuations measured simultane- 


ously are discussed. 





J tallige 


Jill 


aves 


TABLE OF CONTENTS 


PH ESC MOINO coco See eee = eS 8 
Pee A OKIGROUN Dae eee 8 
B. OVERVIEW ------------------------~------------+ Le 
C. PHYSICAL DESCRIPTION ------------------~----~- 13 
D. DESCRIPTION OF SENSORS -----------~-----—---— Pi 
1. General --------------------------------- el 
2. Acoustic System ------------------------- 22 
3. Temperature System ---------------------- 26 
I, AAA WAG UU SIEUSON) MRS SS Se 27 
OUR AWRY STS soc Se Sil a a Ee 31 
fae STR Eee Ri CORDLNG ———————————— ee 3h 
2, 9 DING IE ay NON Sea oe SS ee 31 
8, COWUWERS ION Gae8565555428se6s—.. 2.223565 ==sSSs— 35 
D. PULSE HEIGHT MEASUREMENT ----~-----------~---- 35 
E. SELECTION OF RUNS FOR DETAILED ANALYSIS ---—-— 36 
F, STATISTICAL ANALYSIS ------------------------ 39 
G. SPECTRAL ANALYSIS --------------------------- 42 
H. DISPLAY OF RESULTS ------------~-------------- yy 
RESULTS OF ANALYSIS ------------------—----~------ 45 
A)” TRMPERATURE FLUCTUATIONS --------——--------—- 45 
Bees TlC ws RESURI Oo OF PULSE HEIGHT 
A, POR) Se ee eee 46 
C. SPECTRAL ANALYSIS RESULTS ------------------- 48 
D. DISCUSSION ------------~---~-----~------------ 52 
CRITIQUE OF EXPERIMENT ~--~---------~------------- 62 





ee ee 62 


A 
Pe UNeEN CTH Meee eee 63 
Gee IRECORDLNGMOn = DAGNEeE === ae 64 
DMEEHASE lig RUNSE@e tees 2-=ao-e--2<4—-5------5-- 65 
EMEOREEEECTED SIGNM@c) oe me o=e==-—-—------_ > 66 
MTS CRKERP [NG =o eee 67 
V. SUMMARY, RECOMMENDATIONS, AND CONCLUSIONS ----- 68 
APPENDIX A: PULSE HEIGHT TIME SERIES AND 
ANLYSTS RESULTS —o ee eee 72 
EASE? Gime Tisai Gee eee pl 125 
Mim TD ISTRIBUTMLONULIST -------=seeeees-o--- = 128 
TOG gays 130 





i 


imi 


10 


LIST OF TABLES 


Basic Data on Analyzed Runs -------------~-------- 37 
Sea and Weather Conditions During Each Run ------ 38 
Summary of Statistical Analysis Results --------- 47 


Summary of Selected Results from 
Spectral Analysis ~~----~--~------~~--~~~--~~~--~~-~-~----- Te) 





LIST OF EIGURES 


Map of Auckland-Leigh area in New Zealand -------- 
Map of experiment site near Leigh, New Zealand --- 
Schematic of acoustic propagation equipment ------ 


Typical sample of microscale temperature 
fluctuations and simultaneously observed 
2G OUYCREL GOHAN STS 1G IS Iba gfe) NE la 


Expanded view of acoustic pulse sequence --------- 


Graph of coefficient of variation versus 
number of bad pulses ----------------------------- 


Graph of bad pulses on Hydrophone 1 versus 
those on Hydrophone 2 ---------------------------- 


Histogram from Pulse Height Analyzer for 

pulses recorded during the 5 minute period 

shown on temperature record: low thermal 

activity ----------------- ------- 


Histogram from Pulse Height Analyzer for 

pulses recorded during the 5 minute period 

shown on temperature record: high thermal 
activity ----------------------------------+------ 


16 





ACKNOWLEDGEMENTS 


The author wishes to acknowledge the valuable guidance 
and encouragement of Dr. Warren W. Denner from whom, 
through countless hours of stimulating dialogue, I learned 
the value of innovative and exploratory thinking. The 
assistance and comments of Dr. Edward B. Thornton and Dr. 
Robert H. Bourke are greatfully meee ciated. 

Particular note is due the Department of Physics of the 
University of Auckland Pon the cooperation extended during 
the conduct of the experiment. Special gratitude is offered 
to the staff of the Naval Arctic Research Laboratory, 
especially Mrs. Dorothy Underwood, for the gracious hospi- 
tality and multiple resources made available during the 
author's work at that facility. 

The experimental research discussed in this thesis was 
supported through contracts from the Naval Ordnance Systems 


Command (Code Q3C). 





ee Le ee APE RUMEN 


A. BACKGROUND 

Acoustic fluctuations in the ocean first gained serious 
attention as a result of sonar development and acoustic 
experimentation during and immediately after World War Il. 
Variability in the pressure amplitude of an acoustic signal 
received from a constant source was attributed to the 
presence of moving inhomogeneities which continuously altered 
Ica Cawmepaths Torsenergy TO Travel econ themsourcesto the 
moCceCiVer. 

During this same time frame, the development of quick- 
response thermopiles enabled Urick and Searfoss (1948, 1949) 


to measure the micro-scale thermal structure of the ocean 


near Key West, Florida. Recognizing the prominent role of 


temperature in affecting the acoustic refractive index, they 
proposed that the "thermal patches" evidenced in their 
experiments were the cause of acoustic fluctuations. 

Subsequent investigations — theoretical, laboratory, and 
in situ — are summarized through 1964 by Urick (1967). A 
more detailed theoretical development and analysis can be 
hounaedn Gkwarmyzk €1963)4 

Mtnouph waclearer understanding of the nature of the 
ecean's finem@structure is availabletnow, compared to twenty-— 
Tiye years) aAt@omeune equations often usled to predict ta@coustic 
variability are still those developed by Bergmann (1946) ana 


inc mn commeo> Shue O54). Under the condition that 





r << kaé 


Wiecmecrs re —=— rance irom sSOUrcestTO reicehi ver, 
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acoustic wavenumber, 


0) 
tl 


the "patch size", 


thermal microstructure scale length, i.e., 


the acoustic variability, calculated as the coefficient of 


Variation, V, is given by: 





Ve (get 3 


where u is the RMS refractive index contrast of the 


inhomogeneity. 


At relatively long range, such that r >> Kka* 2 


VamlapiL lity 2S. 2iven Dy: 
qe 2 


iE 
V= ( eel Kae) 


the 


These two equations were derived under the assumption of 


a Gaussian autocorrelation function for the temperature 


fluctuations. The quantities V and wu are calculated as: 


<(P ~<P>)“> Z 
(<p>) 


Cee LiheueniaomVariatlon; eae... V is the 


fractional standard deviation of the pressure 


AMO abe Dwele, Ie 





—— — = pt 
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u = e) » Where c= sound speed 
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RMS variation of the refractive index. 


The "patch size", a , is determined from the space auto- 
correlation function of the temperature fluctuations as 
measured by a sensor moving along a line through the tempera- 
ture field. Invoking the hypothesis of "frozen turbulence," 
and allowing the temperature field to be advected past a 
fixed sensor,- "patch size" can also be determined from the 
temporal autocorrelation function of the temperature fluctu- 
ations, (tT), where t is the lag time, provided the advection 
velocity is known. 

Mintzer (1954) has shown that the autocorrelation func- 
tion for the acoustic fluctuations, R(t), should be related 
onthe =auvecerrelation functiongfor the temperature, fluctu— 


ations by the equation 


ditje> KCe) 


Accepting Mintzer's hypothesis, the temporal autocorrelation 
function of the acoustic fluctuations should be identically 
equal to that of the temperature fluctuations. One recent 
experiment (Campanella and Favret 1969) which attempted to 
verify this relationship found the above expression valid 
Hider cared. COnNLvroOlled laboratory conditions. 

Knowledge of the microstructure of the ocean has expanded 


Seti caiulyvetn the pest ten years. While the classical 
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concepts of thermal patches — spherical, spheroidal, or’ 
lenticular — are still valid under some circumstances, they 
have been largely overshadowed by the discovery of a layered 
microstructure. Evidence of this type of structure has been 
observed repeatedly in widely separated areas of the oceans 
(Stommel and Federov 1967; Neshyba, Neal, and Denner 1971; 
Woods and Wiley 19723; Neal and Neshyba 1973; Gregg, Cox, and 
Hacker 1973, for example). However, the degree of structuring 
and lateral extent of the laminae have been found to vary 
considerably and unpredictably. Thus the tenant of local 
isotropy and/or homogeneity is more difficult to support than 
it is for the case of a patchy microstructure. 

In conjunction with the investigation of a multi-layered 
thermal microstructure, several theories have been proposed 
to explain the mechanisms for the formation, maintenance, 
and collapse of the fine structure. Briefly these include 
internal waves (Garrett and Munk 1972), double diffusive 
processes (Stommel and Federov 1967), and billow turbulence 
(Woods and Wiley 1972). 

In addition to the thermal microstructure, whether it be 
layered or patchy, other factors that cause acoustic fluctu- 
ations may be present in the ocean environment. Among these 
are microscopic bubble populations (Medwin 1973), micro- and 
macro—-biologic effects, and both surface and internal wave 
activity (Crease 1963, Barakos 1972). 

While many experiments have been performed to investigate 


acoustic variability, few have simultaneously monitored the 


dea 





ESSeNtwaleveatures Cie Gne environment. Furthermore, fiew, if 
any, measurements of acoustic variability have been made in 
Shallow water. The experiment described in this thesis was 
designed to collect environmental and acoustic data which 
would lead to an analysis of both temperature and acoustic 


é 


variability in a shallow water scenario. 


B. OVERVIEW 

The purpose of the experinent was to obtain simultaneous 
measurements of both environmental data and acoustic data in 
order to gain a better understanding of the nature of the 
environment, to determine the driving forces evoking changes 
in the environment, and to investigate the interaction ~ 
between the environment and acoustic propagation. 

The environment under study was a shallow coastal region 
With a mean water depth of approximately 70 feet. Several 
advantages were immediately gained from the Satetten ont 
this site. The sensors involved could be rigidly fixed in 
known positions on the ocean floor, thus eliminating the 
problem of platform motion which is difficult to accomplish 
in open ocean experiments conducted from research vessels. 
However, the water depth was great enough to offer the oppor- 
tunity to observe, as the Seasons passed, both a well-mixed 
(and therefore homogeneous and near-isotropic) environment, 
atime Li'd cimmele Gmc rivaitrOnmeiew she extent of ECltTher Situation 
was determined by local wind, wave, and swell conditions. 


Finally, the experiment was maintained for a period of 


2 





months in the same location; this allowed the accumulation 
of sufficient data to form a long term time series. 

During the course of the experiment, the following 
parameters were to be observed and recorded: 


1. high frequency acoustic amplitude fluctuations in 
two Simultaneous orthogonal directions 
temperature microstructure data 


mean water temperature 


surface wave heights 


velocity in two horizontal directions 


vertical velocity 


ambient noise 


‘Surface wind 
stage of the tide 


sky and cloud conditions 


Oo CON NW FW PD 


HoH 
HO 


pertinent phenomena, such as Sunrise, sunset, 
precipitation, unusual occurrences, etc. 


te PHYSICAL DESCRIPTION 

The site selected for the experiment was the University 
of Auckland Marine Station at Leigh, New Zealand. Dr. 
Denner of the Naval Postgraduate School was on sabbatical 
leave with the Department of Physics of the University of 
Auckland. The Physics Department at the University had been 
active for several years in the area of acoustic fluctuations 
ieone OCeadiews Dr. Erick Safar had contributed several 
original papers in the area of sound amplitude fluctuations 
(Sagar 1955, Sagar 1957, Sagar 1960). Furthermore, the 
University Wad previously conducted acoustic experiments at 


Leigh (Brownlee 1969); hence, some of the necessary 


lags: 





facilities were already installed, namely, several conducting 
cables and underwater instrument platforms. The Marine 
Station provided the required logistic support including 
laboratory space, boats, shops, personnel, etc. 

The near-shore area where the range was established was 
being set aside as a marine preserve and intrusions on the 
range by fishing and recreational boats were eliminated. 

The area was very biologically active, especially in plankton, 
nekton, and algae populations (Morton and Chapman 1968, 
Taylor 1972). 

Geographically, the Leigh Marine Laboratory is located 
on the mainland across from Goat Island in the vicinity of 
Leigh, approximately 40 miles north of Auckland on the east 
coast of New Zealand (Figure 1). The range was located off- 
shore from the laboratory between the mainland and Goat 
Island (Figure 2). 

The laboratory had served as a center of marine research 
for a number of years and some of the environmental charac- 
teristics of the area were known. In this regard, local 
climatology was an important factor because the experiment 
was planned to span a period of one year. Statistically, 
the diurnal air temperature range was about 6°C in all seasons 
Gi syne) Yedimmcniserange Wasi only slightly less than the mean 
annual temperature change of about 10°C. The warmest month 
was February with an average daily temperature range of 
20 to 24°C; the coldest month was August with a mean daily 


nanigcuOteMUmGemil> CQ oea SUrtace temperatures historically 
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Figure 1. Map of Auckland-Leigh area in New Zealand. 
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Veriledme«rom an lowmon about 23.5°Cein Aucusteto a high of 
about 20.4°C in February. In the winter the water column 
was isothermal to the bottom. Warming of the water column 
as summer approached resulted in slight thermal stratifica- 
tion which reached a maximum in late October or November. 
Due to the shallow depth at the range, temperature difference 
between the surface and bottom waters was not expected to 
exceed 2°C. No permanent ChenmiOominiom ox s tod Bandeasead sonas 
thermoclines, if present at all, were weak. Surface salinity 
values ranged from 34.20 to 35.86% .witha slight increase 
with depth on the order of 0.1 to 0.2% between ene surface 
and 70 feet. Salinities tended to be Slightly higher in the 
winter and spring months, although this trend was not well- 
established. Light to moderate winds were to be anticipated 
for the majority of the time with speeds seldom exceeding 
18 knots. Currents in the immediate vicinity of the experi- 
ment site were irregular in direction and speed, the primary 
contribution arising from tidal influence. Mean tidal range 
varied from 5 feet during neap tides to 8 feet during spring 
Lldese 

Geologically, the coasts currounding the site were 
primarily greywacke and sedimentary type rocks. The rugged 
coastal topography extended approximately 100 to 150 yards 
offshore where it changed abruptly to a fairly level sandy 
DOfeeoMm UP@iewowehsthe instruments were situated. The bottom 
sediment was friable fine sand of mixed siliceous and 


calcareous origin. High surface wave activity had the effect 


ays 





of agitating the bottom sediments and distributing fine 
particulate matter throughout the water column in a short 
period of time. 

The laboratory was the location of the majority of the 
shore-based equipment. Twenty-one conducting cables from 
the laboratory terminated at a secondary shelter, referred to 
as the "Caboose," where they were joined to three seven- 
conductor armored cables. The three armored cables were laid 
over the coastal cliff and out to approximately 250 yards 
offshore where a set of underwater junctions was located. 
This site became known as the Bird Cage Site, an identity 
referring to the metal protective cages in which the armored 
cables were each joined to a four-conductor cable. Cables 
of the latter type comprised the remainder of the range 
wiring. One of the two hydrophones used in the acoustic 
propagation portion of the experiment was placed at this 
site atop a 20 foot mast; water depth was about 40 feet. 

About ®3l2 yards seaward from the Bird Cage Site was the 
micea Site. <A 20 foot high metal tripod had been erected 
at this "site as an instrument platform for a previous 
experiment. The tripod became the main sensor platform for 
this experiment. It was situated in approximately 80 feet 
Glawaver vice anmanitiial helghnt ef 20 feet@above thesbottom. 
Ae 20) foot@eexvenstion of PVC (Poly Vinyl Chloride) pipe was 
added to the tripod to elevate the sensors to the center of 
The water column. On the tripod mast were mounted the 


acoustic Projector, tvemperature and temperature microstructure 
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sensors, and ambient noise hydrophone. Approximately 100 
feet seaward of the tripod, one Aanderaa current meter was 
moored 20 feet above the sea floor. 

tmecad#ron LO tne hnydrephone placed at “the Bird Cake, 

a second receiving site was established at the base of an 
underwater reef extending northeast from Goat Island 
(Figure 2). Water depth at this location — the Reefesite — 
was approximately 65 feet; the hydrophone was mounted atop 
a 40 foot mast secured to the sea floor. 

By positioning the hydrophones in the manner described, 
two near-orthogonal acoustic paths of approximately the same 
length (actual path difference, determined acoustically, was 
30 meters) were established. The Tripod-to-Bird Cage path 
was oriented in the same direction asmene predominant swell, 
mee, .rOMn Ghewmersheast to ene southwest. Thus, the effects 
of simultaneous acoustic propagation across-swell' and along- 
Swell could be observed. Additionally, if the water were 
Stratified, the effect of internal waves could also be 
observed if present. 

A second Aanderaa current meter was moored 20 feet above 
the bottom at a location approximately 100 feet east of the 
Reef Site. The wave and tide sensor was situated in an 
exposed ameemabout 110 fect south of the Tripod Site. 

MING ti@eniecweducpoment and cablies Had to be anstalled by 
G@rGgers Overarame irom a le foou skitf, In addition to the 
pre-existing cabling mentioned above, new cables had to be 


Mistalled rome cnc bird Care Site to the Tripod, to the Reef 
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Site, and to the wave and tide sensor. These feeder cables 
were terminated in either smaller, lighter cables that could 
casi ebe proughtueco themsurface for sensor connection, or 
in underwater housings in which the instrument connection 
was made. 

Several significant problems were encountered in estab- 
lishing the range. The small boat available required that 
all of the equipment be relatively light. The use of divers 
to install the equipment further dictated that size and 
weight be limited to manageable dimensions. Virtually calm 
weather was mandatory to insure a satisfactory installation. 
Many of the sensors, particularly the thermistors, were very 
fragile and required extreme care in handling. 

Fouling and corrosion proved to be constant problems. 
Experience showed that only a few days! data could be 
Obtained from the thermistors before they had to be cleaned 
by a diver. Biological fouling grew at a rate of several 
millimeters per week on exposed unprotected areas. Stainless 
steel bolts available in New Zealand corroded through in only 
a few weeks. 

imewoceurrence of large waves, either locally generated 
Sue owel | ironmmamdiistant Storm, Occasionally resulted in 
damage to the underwater equipment. Such damage mainly 
eOnsisted of cable displacement and consequent stress or 
eines meeiom MCh LOn  corminals, thus necessitating 


inspection and repair by divers. 
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D. DESCRIPTION OF SENSORS 
ioe CC TIe Tica 
At the-commencement of the experiment in September 
1972, the following sensors had been installed: 
Teeood, oLte: 


. mean temperature thermistor 
Lempera vurce micros vrucuures Therm StOrsarray 


» acoustic projector 


1 

2 

3 

4, ambient noise hydrophone 

5. three components of ducted current meters 
6 


. Aanderaa current meter moored about 100 feet 
DOm@ioe S gO lair 1 ped 


7. pressure sensitive wave and tide sensor located 
about 10 @elet South of theltripod 


Reef Site: 


1. hydrophone 


2. Aanderaa current meter moored about 100 feet east of 
the hydrophone mast 


Bird Cage Site: 
1. hydrophone 

Tne wnorizonvtal components er the ducted current meters 
failed within the first few days of operation. The remaining 
vertical component operated successfully for the duration of 
the experiment. 

All sensors, with the exception of the Aanderaa 
current meters, were hard-wired to the shore station for data 
recording. The Aanderaa meters, which were Savonius rotor 
type meters, digitally recorded current speed and direction 


Tiverially em@emagsnicuae tape for later retrievals. 
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In addition to the underwater measurements, meteoro- 
logical observations were made by meesonnell at the Marine 
Laboratory under a program separate from this experiment. 
Standard atmospheric parameters, including air temperature, 
dew point temperature, barometric pressure, wind velocity, 
precipitation, and insolation, were monitored. 

Two unique aspects of the sensor suite were the 
acoustic and RET Ont ihe components. Early in the preparation 
for the experiment, two special problems were recognized. 
The first was the requirement to measure temperature fluctua- 
tions to an accuracy of at least * 0.01°C; the second was the 
recording for analysis purposes of short ultrasonic pulses 
received at two sites. Temperature was a problem because of 
the sensitivity required over a large temperature range. The 
acoustic pulses were a problem for two reasons: 1) the tape - 
ifeeorder could motu record directly thejhighn frequencies used 
without excessive tape usage, and 2) direct recording of the 
received pulses would require an extremely fast digitizing 
system if the data were to be analyzed on a digital computer. 

systems were designed to cope with these problems and 
are described in the following Sections. 

ea. AcCOUSstIC System 

Mibpwocte diagram of tGhegequipment involved in the 
acoustic Propagation portion of the experiment isjshown in 
Figure 3. Pwo Wavetek oscillators were used to initiate an 
acoustic pulse at the desired frequency. The first oscillator, 


which established the pulse repetition rate, triggered the 
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second oscillator which produced the pulse. The Second 
oscillator, in turn, triggered the power amplifier and 
Simultaneously input the trigger signal to a specially 
designed sample and hold circuit. The power amplifier drove 
a horizontally omnidirectional projector at one of the two 
Sementat frequencies, 35 kHz or 65 kHz. The signals 
from the hydrophones were passed via cabling back to the 
shore station where they were input to the dual channel 
sample and hold circuit. 

This circuit, designed and constructed by personnel 
of the Physics Department at the University of Auckland in 
New Zealand (Ash 1972) operated as follows. The trigger 
Sslgnal from the second oscillator initiated a sample inhibit 
instruction in the sample and hold circuit to prevent the 
circuit from being falsely triggered while the acoustic pulse 
was traveling down range. The pulse output from the hydro- 
phone preamplifier, approximately 2 volts peak-to-peak, was 
amplified in the circuit to approximately 8 volts peak-to- 
peak. By the time the pulse arrived at the circuit, the 
Sample inhibit signal had been cancelled. The first cycle 
from the arriving pulse that exceeded a pre-set comparator 
level enabled a 200 msec sample hold component. At this 
point the peak-to-peak voltage value of the pulse cycle was 
held for 200 msec. The 200 msec output signal was then 
applied to a tape recorder. Approximately five cycles of 
the pulse were required to initiate the sample and hold 


process. After a pulse was identified, the sampling circult 


al, 





was disabled by a logic Sub-circuit to avoid triggering of 
the system by noise or reflections between direct path pulse 
arrivals. 

Pulse durations of 1.0, 1.5, and 2.0 msSec were used 
in various phases of the experiment. However, this did not 
affect the operation of the sample and hold circuit since 
pulse height was determined in the first five cycles of the 
pulse train and (eter arrivals were ignored. 

Pulse repetition rate as determined by the first 
oscillator also varied, with nominal rates of 1 and 3 pulses 
per second in use at different times. Problems associated 
with the higher rate will be discussed in Chapter IV. 

-The three transducers used in the experiment were 
basically identical. Each transducer was comprised of a low 
frequency and high frequency section, resonant to 35 kHz and 
65 kHz respectively. The transducers were Model ITC601D 
manufactured by the International Transducer Corporation of 
Goleta, California. 

Peripheral equipment consisted of an oscilloscope 
which was used to monitor the performance of the sample and 
Mola Girculvumana TO ObServe=the degree of fluctuation 
Securrineg suring wa run. sla thes latter stages of the experi- 
ment, a pulse height analyzer was also placed in line to 
obtain immediate analysis of height distributions. 

Range geometry and water depths were such that 
surface-reflected and bottom-reflected signals arrived at the 


hydrophones sufficiently late enough to insure that only 





direct path measurements were made. During the initial 


phase of the experiment, however, occasional interference 


from surface-reflected arrivals did occur at the Bird Cage 


Seite under conditions of low tide and an acoustic frequency 
of 35 kHz. A satisfactory solution to this problem was 
Obtained” by Towering the mast at the site by 5 feet prior to 
the commencement of later phases. 

3. Temperature System 

Glass encapsulated, fast response thermistors were 
used for measurement of gross temperature and température 
microstructure. For the gross temperature thermistor, the 
Wheatstone bridge components were chosen so that the bridge 
balanced at 14.5°C, the center of the anticipated temperature 
range of 12°C to 17°C. A temperature range of this magnitude 
was expected due to the long period of time planned for the 
conduct of the experiment. 

With the high resolution required for the micro- 
structure measurements, a small change in gross temperature 
would have quickly moved the microstructure bridge off 
balance if a conventional Wheatstone bridge arrangement were 
used. Thus, high resolution would have been available only 
over a very limited temperature range. This problem was 
overcome by using two equal resistors in the ratio arms of 
the bridge and two matched thermistors in the opposite arms. 
tiene leimcrtiasadded tO one sol Whe thermistors ensured that 
only slow temperature changes were sensed. The other therm- 


iacorm imine ad) | aceio arm Sensed both ~ross temperature and 
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temperature microstructure. Thus, the bridge output was 
proportional to microscale temperature fluctuations but had 
the ability to stay balanced over a wide range of gross 
temperature. The stock of thermistors available was 
calibrated in a bath from 11°C to 25°C and proved to be very 
well matched and therefore suitable for this type of 
arrangement. 

For the microstructure measurements, a temperature 
change of 0.1°C was set equivalent to a 1.0 volt change in 
the bridge output; the micros ii@ccare thermistors were 
essentially linear over a range of + 0.1°C from the zero 
point. The gross temperature range of 12°C to 17°C was 
equivalent to a voltage output range of -2.0 volts to +2.0 
werpsewith OsGuvonwseree.Lstered at 145°C. 

To make efficient use of the limited conducting 
cables available, temperature data signals were frequency 
multiplexed in the’ underwater equipment and then transmitted 
to the shore equipment in the laboratory where the signal was 
demultiplexed by filters. The carrier frequencies for the 
two signals were spaced sufficiently far apart to insure 
that rejection of the unwanted signal by the filters was 


effective. 


BE. DATA ACQUISITION 
Preparation for the experiment began during the summer 
eo ln! s sidiciucdced renovation ot the Tripod Site and 


underwater surveys to determine the best locations for 
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additional instrument masts and optimum routing to bring 

the armored cables ashore. Placement of cables, masts, and 
- instruments was executed at various times through the summer 
Of 49/2. 

Phase I of the experiment with all sensors operating was 
conducted during the period 26 September to 19 October 1972. 
Following a six week break when alterations to the mast at 
the Bird Cage Site were made, Snese Il was executed in the 
first two weeks of December 1972. Phase III of the experiment 
took place during the first week Cte Judy o7s., Thesruns 
analyzed in this thesis were all made during Phase Il. 

In order to establish a comprehensive data base, it was 
planned to take measurements on an hourly basis and whenever 
Significant phenomena were observed. During Phase I, this 
schedule was put into effect. However, data storage require- 
ments dictated that the original schedule be modified to 
acquisition of acoustic and environmental data four times 
daily and during unusual occurrences (e.g., high swell, high 
winds, abnormally calm sea, etc.). The four periods specified 
were sunrise, midday, sunset, and midnight. This schedule 
was maintained during the latter stages of Phase I and during 
Phase II. The schedule for Phase III was further modified 
to provide for as near-continuous data collection as possible 
during Ghewewor days that spanned thisspnase. 

Because oepne expansive data accumulation anvacipacted, 
Pinelenaoasaurime, Phases I and fl was planned to be nominally 


LOstciniiLesmmeioewas calculated that @ run of this relatively 





short length would be long enough to satisfy the requirements 
for statistical and spectral analysis of the fluctuations in 
- the temperature and acoustic data. On occasions when extreme 
variability was observed, run length was extended to twenty 
minutes. During Phase III, run length was arbitrary, with 
the shortest run being 27 minutes and the longest 77 minutes. 

In order to investigate the frequency dependence of the 
acoustic fluctuations, runs were generally conducted in 
series of two, alternating the high and low frequencies in 
each series. Due to the mechanical manipulations required 
to change frequency, the two runs in a seriesS were Separated 
by a period of about 8 to 20 minutes. 

Two magnetic tape recorders were in use during Phases I 
and Il of the experiment. One was a four channel Hewlett-— 
Packard Model 3960, with a tape capacity of 1800 feet of 
1.5 mil tape; the second was a seven channel Ampex FR1300 
model with a capacity of 2500 feet of 1.5 mil tape. On the 
Ampex tapes mean temperature, temperature microstructure, 
pulse heights from both hydrophones, and the current data 
from the three ducted current meters were recorded. Wave 
height data, turbulent velocity data from the vertical compo- 
nent ducted current meter, multiplexed temperature data, and 
pulse height data on a time-sharing basis were recorded on 
the Hewlett-Packard re covdecrmue.ne firsguehna li sor ga run was 
dedicated to Hydrophone 1 at the Bird Cage Site, and the 


second half to Hydrophone 2 at the Reef Site. 
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Aanderaa current meter Tapes wereyprocessed after 
retrieval by displaying the data on a calibrated strip 
chart. Except for a few runs in Phase II, ambient noise 
was not recorded, but was monitored during each run. Also 
in use was a four channel strip chart recorder which dis- 
played wave and tide information, gross temperature, and 
temperature microstructure data on a real time basis. 

Log-keeping involved recording of dates, times, tape 
recorder channel allocations, significant wave height (Hy 73) > 
Stage of the tide, and a summary of weather conditions. 

Wave height and tide information were determined to the 
nearesi:,0.. feetefrom the’ strip chart. 

Because the experiment was a joint venture between 
personnel from the University of Auckland and from the Naval 
Postgraduate School, copying of data ‘tapes was required 
after each phase to make the data available to both parties. 
Selected tapes from the Ampex and Hewlett-Packard recordings 
were transcribed onto one inch, fourteen track magnetic 
tapes with the use of a Sangamo 3562 tape recorder. Addi- 
muna lly ecopsles of LThemexperimeng, logs ana of somewot the 
current and wave strip charts were returned to the United 
otates for analysis. 

During Phase III, the Sangamo recorder was used for 
initial data acquisition and dubbed tapes were provided to 


Ene Universroy of nuckland: 





IT. DATA ANALYSIS 


A. STRIP CHART RECORDING 

The first step in the analysis procedure was to survey 
the transeribed data to determine qualitatively the nature 
SL the experimental runs available. Initially, this) was 
attempted by viewing the data on a dual-trace oscilloscope. 
However, only two channels could be viewed simultaneously 
and no permanent visible record was produced. 

The magnetic data tapes were then played back on the 
Sangamo tape recorder with the output displayed on an eight 
channel Brush strip-chart recorder. The two tape channels 
containing hydrophone information were input directly to the 
Brush recorder. Each channel containing temperature data 
was first low pass filtered through a Krohn-Hite Model 3340 
filter at 20 Hz to eliminate high frequency noise that may 
have contaminated the signal. The filtered data were then 
displayed on the Brush recorder chart. 

A typical section of pulse height data and temperature 
microstructure data is shown in Figure 4(a). Figure 4(b) 


shows an expanded view of a typical pulse Sequence. 


Bee DIGITIZATION 

After screening the data for noise and run length, thirty- 
twee runsmmade during thesperiod 30 Septemberm@through 13 
eroocr mmo femmca@asliVe ruse ton thesperiod’ 2 through 3 July 


1973 were selected to be digitized. None of the runs from 
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Phase II of the experiment were chosen. Four channels of 


acoustic and temperature data were digitized. Digitization 


was accomplished through the use of a COMCOR 5000 analog 


coneucer. 7a sclentific Data Systems XDS 9300 digital computer, 
and a Naval Postgraduate School analog-to-digital conversion 
program. All of the above were available in the Electrical 
Engineering Computer Laboratory of the Naval Postgraduate 
senool. 

ingthe digitization process, the most critical factor to 
be specified was the sampling rate. Based on a requirement 
to accurately define the pulse height on the digitized 
record, a sampling rate equivalent to 256 samples per 
channel per digital record was selected. This rate translated 
to a sampling time increment of 0.04 seconds between samples. 
Since the length of a pulse was 200 msec (0.2 Seconds), the 
Sampling rate insured a minimum of five samples on a normal 
Pulse height and a maximum of six samples if the first sample 
on the pulse happened to occur exactly at the leading edge. 
Each hydrophone channel was input directly from the Sangamo 
tape recorder to the analog patchboard where the signal was 
amplified 25 times. 

The Nyquist frequency associated with the sampling rate 
Wasmgaven by Che) equation, 


= 41 
i ~eN 


where At = sampling time increment. 








For the increment chosen, 0.04 sec, the Nyquist frequency 
was 12.5 Hz. This frequency was above the highest frequency 
of significant fluctuations expected in the temperature 
microstructure data and was high enough to minimize aliasing. 

Temperature data were played into the input of the analog 
patchboardg where a plus 10 gain was applied. The amplified 
Signal was then low pass filtered at 50 Hz se © iminate 
locally generated 60 Hz noise that was apparent in initial 
trial runs. Because tape playback during digitization was 
conducted at 16 times real time, the actual temperature data 
was, in effect, filtered at 3.125 Hz. The filtered signal 
was then input to the patchboard where a second plus 10 gain 
was applied. Total amplification of the temperature data 
was 100 times. Of several methods evaluated for filtering 
and amplifying the temperature data, the above method 
appeared to be the most practical and efficient in eliminating 
the locally introduced 60 Hz signal. 

In addition to the above procedures, biasing of the 
amplified signal was required to keep the signal voltage 
within the = 100 volt range of the analog computer. The 
COMCOR 5000 provided bias adjust on each channel. 

Bach anealos run was converted to a digital file of 
Giscrete samples. These samples were written on Seven track 
magnetic owe i Mmeceitet base movdlc On we Lachaiile containdd 
PT maoMimechanielemet digiti zea data appiitcable to that run. 
Data in a file were written in blocks of 1024 samples called 
digital records; each record spanned 10.24 seconds of real 
time Pate 
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C. CONVERSION 

Continued analysis was to be performed on the IBM 360/67 
 Gigital computer at the Naval Postgraduate School. However, 
it was first necesary to make the digitized data compatible 
with the IBM 360 system. The conversion process involved 
changing the notation from octal base to hexadecimal base, 
and rewriting the data on a nine track magnetic tape. The 
conversion was made using programs available at the Naval 
Postgraduate School computer center and described in a 


current technical note (Raney 1973). 


Pee cUloe HE TCRTeExATRACTION 

Of the 256 samples per channel in one digital record 
which represented 10.24 seconds of data, only about 50 of 
the acoustic channel samples actually fell on a pulse height 
peak. Furthermore, these 50 samples accounted for only 10 
Pulses. In order to conserve core storage, it was necessary 
to compute and store only the pulse height value, rather than 
read into memory an entire file, consisting of an average 
58 digital records. 

MW program was written to read from the tape’ only the two 
channels containing pulse data. When one digital record had 
been read in, the pulse heights for each acoustic data channel 
were computed and stored. Upon completion, the next suc- 
ceeding record was read from the tape and processed. This 
sequence was continued until the entire file had been 


analyzed. 





Programming precautions were taken to insure that no 
pulse heights were missed on either channel and that only 
actual pulse heights were retained; i.e., spurious noise 


spikes were omitted. 


E. SELECTION OF RUNS FOR DETAILED ANALYSIS 

Of the 35 runs in Phase I that were digitized, three were 
unacceptable due to excessive random noise that was apparently 
acquired during digitizing. From the remaining 32 runs, 
CWelves Were selected for detailed analysis. The selection 
criteria were to form a suite of runs that were representative 
Of all typical Environmental conditions, and to choose runs 
that for one reason or another (e.g., frequency, time of day, 
time sequence) were suitable for comparison with each other. 

No runs from Phase III were considered due to low apparent 
Varian ty=and difficulty with the? data resultine from 
unexplained interchannel feedback somewhere in the data 
acquisition stage. 

Pertinent details of the twelve analysis files are shown 
rime! abates | ebecause ef digitizing requirements, thertfile 
liength was always slightly shorter than the actual experiment 
run length. Table II Summarizes the environmental conditions 
present during each run. The comments section of Table II 
was synopsized from the original logs kept during the 


experiment. 


36 





"UOTSSNOSTD 
quenbsesans Uf 04 PerIdeyed se *Sated UF pseaonpuod 
sund JO 39S 944 34n4F9SU090 JTT pue OTT pue *T9 pues Og 





“eG pue 9G Soh PUe Gh SUNY “SUNY peZzATeUY UO e4eq OTseq “IT sTaQey, 
18S SyE°?cG We PASTA EIE = 1E G9 WUE 1 SE 390 ET nial 
L8G Scc’c WOL 89Lc20°T Gt be 0 iaece | 2G) SIL O9LL 
rie l Sgc°S Wé6T gETC6°0O G9 9€02-9T0e 490 ¢ Su 
TT9 SgQ°Te wW6 61026°0 G9 g2C L-O teu 490 ¢ Cis 
SL9 Scl°cG wWe6 G9096°0 G9 Sco Usbiroo 49O L US 
QT9 sgol°cG we6 TL096°0 Ge OT60-0060 4°90 I O09 
629 S9O°€ wWOT 62096°0 G9 9¢c90-9L90 420 I gS 
LT9 SQL°TS wWe6 89096°0 Gt eGG0-cr7S0 490 TL 9¢ 
8 L9 Syrl°cS wW6 89096°0 G9 9n9T-9E9T das ot oy 
ETO SHQ°TT > weé OLT06°0O Ge 8cST-STST des o€ q 
rL9 SQ6°TT wé6 97006°0 G9 DegkeOcce dss ot T 
nL9 -SO6°TT Wo €£006°0 Gt O ea O Cm des 0o€ Gi 
(998 * UTW) (99S) 497 

sos—nd U4SsuUsT VUSUSIOUT (ZH) potdsg (C)icse TIQuit 

JO OT Ta SU TL Adusenbsd gy SwTL 342d uNnyY 
ISQuUNN ‘doy ostndg TeooT 


N 


37 





UuNnYyY UoeY BUTAING SUOTATPUOD ZSUIeSM DUS eBdaS 


“00SO +e Peqtodat ,wyed,, worl suTseatouT AT Tenutquoo 

uaeq peu SpuTM pUe Sep TTS uns ou ussgq pey azeyd S)TT pue OTT sund 

Jo otep ou 6.18Q0990 ET UO “OTT UNI 04 JoTad pepauooar sem uoTJeqTdtosid 
ON ‘*4SamuyTOU 09 4SomM WOTT ATUPeU ,YSaTy,, OF Wed, WOT SUTATCA 

S19M SPUTM “JTLT PUe OTT stmt BSuTpsdaid sumoy gy doy 4SeouwsAo useq pey ANS 





"IT eTGesL 


:o40Uu $,z0uqny 


‘PpUTM Ysety “4SeotaAg g°T 9s a 
‘uTer ATTenbs ‘4semM WOT] PUTM USety 04 aqerepow SyTep ATTesu *4seo7eAC Oa L 0°?¢- neu 
“PUTM 4USTT Adena “HIeEq 9°0 Oca € 
"4Sem WOT pUTM AUSTT Area *(QOBT 98) UMop Suto’ 4snf ung G*0 O° T+ ee 
‘uTeo vas SU4noOSs WOTy JUBTT AJOA PUTM S2SA0CD PNOTO OT/Z 9°0 Q°T- T9 
*TIoems ATaeyugnos ST a4e94s eas *4UBTT ATOA PUTM SISAO0 PNOTO OT/T Powe C7 09 
[*podt2z4 syuq ye SUTATP pantesqo 
O18M SOTeUM OM4 *()+190-1£90) Unt Buppesdons BsuTang] *suoydoapAy astou 
quetqwe uo presy STeUuUM JETTY Sdn Burwoo qsnf uns ‘spnoto ou *purm oN 9°00 Orc QG 
*purmM ou *zTeaTo ANS 6°0 11S ie 9S 
“wTeO ST podt.14 Teou ess SpuypM ATAaqseo 4QUBTT Aaaa Suns 4usTIQ SteatTo AysS G6 0°O 64 
"yqnos 9844 WOTJ PUTM 4USTT ATeA Spamosqoun uns 4nq “ApnoTo ATI1eg 9°0 Q*T+ Qh 
“E¢eT 9e upese 
qno yoreq ums Sunt Jo 4.7e9S 4e pnol[oO Aq pamosqo uns *2daA00 PNOTO OT/E 9°0 Q° E+ Of 
“puypm 4UusTT Adea *SAuUNS ATYSOU SASAOCD PNOTO OT/Z 9°0 G°o+ Gh 
SJZUSUUIOD DUS AdYUAedM C37) e/T 3, Gane) TOQuUNN 


JUSTOH SACM ODTL uny 


38 





F, STATISTICAL ANALYSIS 

In order to perform the planned spectral analysis, an 
accurate determination of the pulse repetition period was 
required. To obtain this, the average time increment, At, 
between pulses (Table I) was computed. This increment was 
the measure of the frequency with which an acoustic pulse 
- "sampled" the medium. 

The time increment, At, was found by computing the total 
time elapsed between the first sample on the first pulse in 
Poul wwe che first "sample on the last pulse’gjin the same 
He est cmeelapyed=time, referred to as@file length an 
Table I, was then divided by one less than the total number 
of pulse neauetnee y 

The analog-to-digital sampling period of 0.04 seconds 
was the maximum error possible between the actual analog 
five #enestcheanmd the digital fille length @s computed above. 
Thus, for a file having a nominal 601 pulses (see Table I), 
the maximum possible error in the computed average pulse time 


increment, Aji, was on the order of: 


Ca eee = = eo = Bot so) seconds 
error 


Am error Gf this magnitude was consideredseaeceptable for the 
precision required to perform a valid spectral analysis of 


maemo L Semel ences . 
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Computation of certain statistical parameters was the 
next step in the analysis scheme. A pre-compiled subroutine 
available in the Naval Postgraduate School IBM 360 library © 
wae UsedesbO Compute the coeffictent of variation, the 
minimum and maximum values and the range of the distribution. 
The pulse height histogram was also computed and plotted. 

It was immediately apparent from the histograms that a 
Small percentage of pulses had either extremely high values 
or low values. However, this result was expected from the 
Brush recordings which showed a random occurrence .of 
"drop-outs" and "spikes". The origin of these anomalies was 
not exactly known, but it was apparent that they did not 
represent real acoustic cienuntieee To eliminate these 
a program*was wriltten to filter the pulSe heights through a 
gate whose width could be varied. The acceptable gate width 
was determined by a trial and error process. The ‘method 
employed in this smoothing process was to compare a pulse 
height with a running mean computed from the previous ten 
pulses. The first ten pulses in a file were tested against 
the mean of the first twenty heights to initiate the’ process. 
The gate was centered around the mean. A bad value was 
Se Mima tedepyasubSsSlClLbubImg an a@rvilicial pulse whose height 
was determined from a linear interpolation between the 
Prececolie oooampulsesand the next Succeeding good pulse. 

PU sceehecigiveswtiitered in each computer run were compared 
with the Brush chart to determine if any food data had been 


filtered or if any anomalies had been missed. This method 


NO 





Gid not eliminate all bad pulses and on some runs one or 

two erroneous pulses escaped correction. Since these few 
represented a very small percentage of the total sample — 
less than 0.5% — no further attempts were made to smooth 
them. Thus the temporal fluctuations in the smoothed 

EL SIC DUG Lom Were. Wath raresexception, the result of actual 
acoustic variability. 

The smoothed data were then-.re-analyzed as discussed 
above. Of all the basic statistics computed, the coefficient 
of variation was the most pertinent. Skewness and le bugieepelilts . 
Which were computed along with other statistics, would have 
been valuable for comparison of histograms but were 
invalidated for this purpose because the class interval 
varied among distributions. In computing the histograms, 
the constant faehor was the number of class intervals rather 
than the width of the class interval; this feature ERE Ks 
characteristic of the subroutine. 

To facilitate comparison of results among runs, the 
pulse heights of each run were normalized to the mean for 
that run. Each pulse height sequence, then, had a mean of 
1.0, and each pulse value was expressed as a percentage of 
the mean. This process had no mathematical effect on the 
SOcliimelcneson Vvariatton. Conversion of pulse heights to 
Standard variables’ was also done so that distribution in 


terms of standard deviations was available. 
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G. SPECTRAL ANALYSIS 

Following procedures detailed in Bendat and Piersol 
(1966), preparation of the pulse height time series for 
spectral analysis started with removal of the mean and of 
any linear trend present in the series. 

An estimate of the autocorrelation function, R(t), for 
each time series was computed using the formula: 

N-m 


ROT) = iE Oy) Xa) 
n= 


KJ 


where N = the number of time samples, i.e., the 
number of pulses in each run 
m= the lag number 
t = the lag time = m(At) 


X = the individual pulse height value 


The values for tT were incremented from tT=0 to a TeReuTIG lag 
value of approximately one-tenth of the file length. The 
increment used was the At determined by the method in 
meocl Lone: | el. 

A Parzen lag weighting function was then applied to the 
BU vee oiler eC lol 1 Omni UC i Ons me Order eEOuiisure a betler estimate 
of the spectral energy density values. 

Thnesautcecorrelation function was Fourler transformed to 
arrive at an estimate for the spectral energy density 
function, G(f). The digital programming was derived from 


ioe Inuemradl #exDression: 
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oo 


Gr) -e4e sf “Rit) coseenft drt 
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where R(t) = the autocorrelation function 


tl 


it frequency 


T lag time 


Estimates were computed for narrow frequency bands over the 
ePamcestron Dota te f. where f is the Nyquist frequency 
applicable to each run. 

Themen SSmecorresari1on Lunct Lon Rt? between Hydrophone 
1 and Hydrophone 2 was computed next, using the equation: 
1 


Ay Cl Net 2 ) ie G 


iH} 


where X pulse height value for Hydrophone l 


y pulse height value for Hydrophone 2 
and N, m, and ‘t)were’ as previously defined. 

As with the autocorrelation functions, a Parzen window 
was USeEGdsEOG smooth the cross correlation function prior to 


computation of the cross-spectral energy density function. 


The co-spectral density and quadrature spectral density 


Minions were mused OnCbua lnmehe seoherence. a » and 
phase, Oy merinctions fOr the pulse height time series: 
° oO 
ce... Cy Ceimoe (") 
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coherence between Hydrophone 1 and 
Hydrophone 2 as a function of frequency 
phase between Hydrophone 1 and 
Hydrophone 2 

co-spectral energy density function 
quadrature spectral energy density 
Prime vl Oe 

Spectral energy density function for 
Hydrophone 1 

spectral energy density function for 


Hydrophone 2. 


Hee Di ShLAYe OF RESULTS 


tic 71 oklowine compurer generaved plovs@were®made on the 


CABECOMP plotter: 


1. relative pulse height time series for each hydrophone 


BULCCor rue lal Ole Lumet tons Ei cGas cachenydrophone 


auto-spectral energy density of the acoustic 
fluctuations for each hydrophone 


4}. phase and coherence between Hydrophone 1 and 2. 


Except for one instance (Run 73) where the file length 


was longer than usual, all runs were plotted using identical 


Scobling factors. 
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i neo Uke On ANATYSIS 


A. TEMPERATURE FLUCTUATIONS 

Detailed analysis of the temperature microstructure data 
Was miOmepcrLOrncd gine the course Gf this thresis. However, 
the Brush chart recordings previously described afforded an 
opportunity for a’ qualitative evaluation of the micro- 
structure activity present during each run. Of the twelve 
runs under discussion, four showed temperature fluctuations 
that were considerably more active than the others. The 
remainder of the runs occurred under temperature conditions 
varying from negligible to moderate activity. The most 
active runs were numbers 72, 73, 116, and 117. Representative 
samples of the twelve temperature records are shown in 
Appendix A. 

The temperature microstructure fluctuations weve 
evaluated using the thermistor calibration factors, tape 
recorder gains, and the gain and chart speed settings for 
the Brush recorder. An arbitrary classification scheme was 
employed to assign qualitative descriptors: 


Temperature Fluctuation Classification 


ee 5.6 Negligible 
nOS CML O 9.0.2°C Moderate 
27 eo © High 





B. STATISTICAL RESULTS OF PULSE HEIGHT ANALYSIS 

Table III summarizes the statistical results obtained 
from the analysis. Reference to Tables I and II supplies 
information relating the neonate frequency, time of day, 
and environmental conditions to these results. 

The coefficients of variation ranged from a low of 2.03% 
(Run 56) to a high of 15.53% (Run 116). The five highest 
coefficients occurred in runs of relatively high thermal 
ectivity. ln general, the variation observed on Hydrophone 2 
was Slightly higher than that on Hydrophone 1, the exceptions 
being Runs 45, 46, and 72. However, in only four runs (60, 
73, 116, 117) was the difference between the two hydrophones 
greater than 2%; three of these four runs were also classified 
oo neaving higher microstructure content than the others. 

Significant differences between the high frequency runs 
and low frequency runs were more difficult to determine. For 
runs that were conducted in high and low frequency pairs (see 
Table I), the coefficients of variation for both Hydrophone 
lL and 2 during the high frequency run were generally slightly 
higher than during the low frequency run. However, since 
the two runs in a pair were conducted 8 to 24 minutes apart 
(Table I), the difference may not have been due to frequency 
change alone. 

Table III also has two sections pertaining to the range 
Mime ieeClsurepubtons Of pulse heights. The first of these 
sections lists the extremes of the distribution in terms 


of standard deviation. The second shows the range in terms 
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of percentage deviation from the mean. Information about 
the shape of the distribution can be inferred from these 
statistics. <A relatively large o range coupled with a rela- 
tively small percentage deviation indicates that the dis- 
tribution is closely packed about the mean; this is reflected 
in a low coefficient of variation. The converse situation 
implies a relatively flat, dispersed distribution with a 
consequently higher coefficient of variation. The coeffi- 
cient of variation can be obtained by computing the ratio 

of relative percentage range to standard deviation range. 
nevevcremipne Jcrining  =equaviton given am Chapver I was used 


emcee NremGnc COCtMmaCTLCNUS Pimemen lable rt 


C. SPECTRAL ANALYSIS RESULTS 

Computer generated plots for the twelve analyzed runs 
are presented in Appendix A. Each series of four plots for 
each run contains the time series of pulse heights plotted 
relotave Som Lies Meakeeuhe auvocerrelataon Tunerions for 
each time series, the auto-spectral energy density functions, 
and the coherence and phase functions determined from the 
cross-spectral analysis between Hydrophones 1 and 2. 

Table IV summarizes the lag times determined for three 


ae Ehovsl OOD) 


values of the autocorrelation function (0.5, e 
approximate values of strong periodicities observed in the 
HmecerrelarlommruncuLlons, andavne maximum value of coherence 


attained between the two hydrophones. 


ike 





Autocorrelation Lag Times 


Run Hydro- z (sec) 2 Apparent 
Number phone On Ve 0.0 Periodicity 
4s 5.5 8.1 ha garg None 
2.0 Eee) So aS None 
46 7.4 9.8 28.8 None 
1.2 2.0 19.8 8 sec 
48 ec) 4.0 19.3 None 
TT il 17.6 = None 
hg 13.9 Ios 28.8 None 
5.3 WMlab 29.8 None 
56 0.8 ilgeits aaa None 
Le flier Risse 8 sec 
58 0.9 1.9 19.2 None 
ileus 2 all 11.0 None 
60 4.8 16.3 _ None 
eel loa = 7 sec 
61 0.8 a PO 6 sec 
1 Tago’ 3.5 6 sec 
Te - -~ None 
1226 545 AAS a None 
13 2 AL 6.5 31.3 None 
10.9 1 ie _ None 
16 des a eee Pu oe 4h sec 
2 ee = 7.5 sec 
ia, less) lee ho 4 7.7 sec 
2.8 dite _ 6.5 sec 
Table IV. Summary of Selected Results 
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The autocorrelation functions exhibited a wide variety 
Semenareacveristics Which were related to the acoustic fluc- 
tuations, and consequently to the state of the medium at 
the time the run was conducted. In arun when the temperature 
microstructure activity was low to moderate, evidenced by 
the figures in Appendix A, and the wave activity was low 
(Table’I1I), the two autocorrelation functions were roughly 
comparable. For runs that occurred in pairs (Table I), 
no reliable trend was detected to indicate significant 
difference due to a share’ IM acoustic PREC uenC tim, When 
periodicities were apparent in the autocorrelation functions, 
Hic m~vcre steoOntemeon Hydrophone 2 vnem on Hydrophone 1. 

Runs 7e and 73, conducted two moweee apart during periods of 
relatively high thermal activity, did not evidence any 
sLrong@periodicity, uniike Runs 116 and?1l/, conducted under 
apparently similar thermal conditions. Significant wave 
Neient ae tormtnese LWORMSCLS Gm rule were 0.6 feet and 1.8 
heewmecespectavely, Forjruns®where the® autocorrelation 
function showed a marked periodicity, a distinct peak in 
the spectral energy density was evident. However, many of 
the peaks in the energy spectrum were not manifested in 
meadtily detanablle periodicities in thes autocorrelation 

AUN OM Onmalon at 

A! vit WOmdUmocorre latlonetunctdons exhibited 2a derinite 
trend towards zero with increasing lag time. Hydrophone l 


for Run 72 was decaying very slowly, but had not reached 





a value of 0.5 at the maximum time lag of 55 seconds. In 
; Run 49, Hydrophone 1 decayed continuously from Tt = 0, 
Peotatne the zero axis at 30 seconds and continued to decay 
until it exceeded the scale capability of the plot at a 
value of -0.34 and tT = 43.2 seconds. However, the computer 
printout showed that the autocorrelation function reached 
a minimum of -0.38 at 51 seconds and was increasing towards 
zero at the maximum lag time of 58 seconds. 

For all runs, the coherence between Hydrophones 1 and 2 
was very low with an average value on the order of 0.1. 
The maximum coherence listed in Tabie III occurred on a 
coherence peak that was not representative of the entire 
run. A maximum coherence of 1.0 would have indicated a 
linear relationship between the fluctuations on Hydrophone 
1 and those on Hydrophone 2; conversely, a coherence function 
of 0.0 would have meant that the fluctuations were unrelated. 
POmVeatiimermnealtatem—valves ,  Lnree —e@ocrt bile ssuinvegauerne ero ulalte. 
have existed [Bendat and Piersol 1966]: | 

1. extraneous noise was present in the measurements 


2. the relationship between the fluctuations on each 
hydrophone was non-linear 


Ss. the fluctuations on one hydrophone were not driven 
entirely by the same mechanisms as the fluctuations 

One tne Orner hydrophone 
Vittoria tanmom mice CXDeCriImental Daramemers, a non= 


linear relationship between the fluctuations was the most 


plausible reason for the low coherence. 





The mwieasCmmunceLon=t Or all@yuns 1S #veryeecrracic, a 
Zeiurpvner indrPcation that .a simple linear relationship did 


not exist between the two hydrophones. 


D. DISCUSSION 

The success of the analysis of the acoustic data was 
predicated on the reliability of the sample-and-hold circuit 
described in Section 1.D.2. Therefore, due scrutiny was 
given to the number of bad pulse heights that were detected 
Peeeunes SMOOLhIng sarecesis outleine dn Section 22.m. Jf the 
number of bad pulses were proportional to the acoustic 
variability, the probability that such pulses may have been 
actual acoustic fluctuations is increased. Plotting the 
coefficient of variation versus the number of bad pulses, 
Figure 5, did not support a relationship between the two 
parameters. Data for Figure 5 were taken from Table III 
and all points were plotted except those for Run 73 which 
was twice as long as any other run. 

TnhewauMnber ven Dela pullses’, However, did seem to be related 
to frequency. There were significantly more rejected pulses 
atwo5 KHZethan at 35 kHz This suggested that the sample 
aid HOlamerrcult wasmencounterings trouble@im= ‘Tollowing” 
the 65 kHz wave form. However, Hydrophone 2 showed 2 to 
BeuuMcoe Ome bad pulses than Hydrophone 1. FA plot of the 
number of bad pulses that occurred on Hydrophone 1 against 
the number on Hydrophone 2 for both frequencies, Figure 6, 


Pigdieaved scurorm™® relaurdonsnip between the two parameters. 
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Since the hydrophones were identical and operated well 
throughout thevexperiment, the bad pulses could only have 
been generated in the two channels of the sample and hold 
apparatus. The critical question was: Did the malfunction 
nullify the experiment? Evidence to the contrary was ex- 
plicit in the pulse height histograms which clearly showed 
the bad pulses separated from either side of the main 
grouping. Furthermore, the Brush chart time series supported 
the conclusion that, except for the random occurrence of 

a bad pulse, the sample and hold circuit performed reliably 
on all other pulses. 

During Phase II of the experiment, an on-line Pulse 
Height Analyzer was installed to obtain a real-time analog 
display of pulse height distributions. sSample histograms 
from two runs made under conditions of low and high micro- 
structure activity are shown in Figures 7 and 8. The Pulse 
Height Analyzer was storage limited to the number of pulses 
in a five minute segment of the run. The nacnaseons die five 
minute portion of the temperature microstructure record is 
shown for comparison. The histograms obtained by digital 
methods in this thesis were comparable to the distributions 
accumulated by the Pulse Height Analyzer. This promoted 
additional confidence in the sample and hold electronics. 

Pemnobegeprevyilously, the degree omerhermalf¥activity 
varied considerably during the runs from Phase I that were 
Balyzedwumenuns (2, (3, 116, and 117, all conducted unde? 


eonditions of similar temperature fluctuations, exhibited 
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notable differences in acoustic variability. Referring to 
Tables I and II, Runs 72 and 73 were made under conditions of 
very light wind and a virtually calm sea surface over the 
acoustic range. The runs conducted the previous day (Runs 
56-61) also experienced calm weather. From the logs, the 
conditions on 2 October prior to Run 72 were similar with 
calm winds and sea and a high, thin broken to overcast sky. 
The large acoustic fluctuations observed during Runs 72 
fees) Were Indicarive of the microstructure activity. No 
strong periodicities were apparent in the- autocorrelation 
functions and the spectral energy density functions for both 
hydrophones during each run were similar. 

In contrast, the weather and sea conditions (Table IT) 
were different during Runs 116 and 117. The acoustic, as 
well as temperature, fluctuations were large, similar to the 
neMomOnmc sOCLODCY »e HOWCVEr. Cheraucvocomrelation functions 
exhibited marked periodicities and the spectral energy 
density functions had a distinct peak at about 0.14 Hz. 

nyt noueh thelcomparison 2s made®*between two sets of runs 
conducted eleven days apart, one conclusion appears reason- 
able. The acoustic fluctuations were induced by at least 
two mechanisms in the medium, the temperature inhomogeneities 
anomvunes ci iect OF msurfacle Wave. activity. This tsituation has 
been observed in other experiments as well [Alexander 1972, 
| for example]. Confirmation of the effect of surface waves 
would require an analysis of the wave spectrum, the data for 


which were not available for this study. However, the periods 





seen in the autocorrelation function (Table IV) are readily 
observed in a surface wave spectrum. A sample calculation 
of the water particle velocity at a depth of 20 feet indicated 
that the velocity effects should have been small for the 
wave heights present. However, the presence of waves can 
affect the transmitted acoustic energy in the wave spectrum 
if inhomogeneities are present, viz., temperature and salini- 
ty contrasts, bubbles, biologics. Then the wave induced 
particle motion would serve to move these inhomogeneities 
in and out of the propagation path in a périodic manner. 
The pulse height time series for Runs 116 and 117 also 
presented evidence of long period fluctuations upon which 
the higher frequency excursions were imposed. The oscilla- 
tion was distinct on the records for Hydrophone 2, but 
barely discernible on Hydrophone 1. Run length was too 
Sfocenyo Gaiane thus fluctuation in the spectral analysase 
although its effect was evident in the autocorrelation 
functions for Hydrophone 2. The ‘ada coefficients of 
variation for Hydrophone 2 in Runs 116 and 117 were caused 
by the non-stationarity of the mean associated with this 
mineral On se Lf CThemcocti ficients forgtheise two runs are 
Mmescounved, Ghee vermavlon forethesset of 12 Prurs® analy zed 
then ranged from 2% to 10%. The apparent period, greater 
Chaviwe: ie mMinuLles, and Gs occurrence In two runs conducted 
eighteen minutes apart indicated the possibility of the 
presence of internal waves. Climatology, discussed in 


Becuylon ).¢., showed a trend for thermal stratification to 





occur in October as’ the summer months approach, thus pro- 
viding the requisite density interface for internal waves. 
Other runs, particularly 60, 72 and 73, showed indications 
Mirlonomperiog Llneuwuations, but not as distinctly as Runs 
116 and 117. Analysis of the temperature microstructure 
data would not be conclusive in evaluating the internal 
wave hypothesis because of the short run length compared 
to the period of fluctuation. A similar situation was 
reported by Shoemaker (1971) in the Arctic. 

The definite appearance of this "wave" on Hydrophone 2 
and its simultaneous absence from Hydrophone 1 were suppor- 
tive of the anisotropy of the range. Further evidence of 
a direction dependence was apparent in the time series plots 
of pulse heights and in differences in the autocorrelation 
functions for other runs, notably those runs where a marked 
periodicity occurred (Table IV). In those cases, Hydrophone 
2 always exhibited a stronger periodicity than Hydrophone 
Pew lineall put three rons (45, 46, 48), the total spectral 
emerey, derived from computing the» area under the spectral 
energy density curves, was greater for Hydrophone 2. The 
apparent cause of the difference in acoustic fluctuations 
Pont wmicm uo Orvmoconal paths was the aniluence of surface 
waves. 

In addition to the possible causes discussed above - 
temperature inhomogeneities, internal waves, and turbulent 
water particle velocity - other factors that nay have 


influenced the acoustic proparation were bubbics, sediment 


60 





particulate stirred into the water column, and biological 
activity, not the least of which were the whales observed 


during Run 58. 
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HY.) 6 CRITIQUE OF EXPERIMENT 


A. PARAMETERS 

Data on several additional parameters would have aided 
in the analysis of the acoustic data. The first of these 
was vertical temperature profile information. No profiles 
were available to determine the presence or confirm the 
absence of seasonal or diurnal thermoclines during the runs. 
Such information could only be weakly inferred from clima- 
tology and current and past weather data pertinent to each 
run. A temperature profiling system was designed and built 
but technical problems and lack of conductors precluded its 
operation. The requirement for a remotely controlled system 
Which would operate reliably from a known, fixed location 
under a variety of sea state and temperature regimes was 
not attainable in the time frame of the experiment. The 
other alternative was to take bathythermographs from a 
vessel stationed on the range during each run, but the 
instrumentation was not available. 

Turbulent velocities measured by the ducted current 
meters were limited to the vertical direction after the 
fai LUE Of ether horizontal components in the initial stage 
of the experiment. Although the remaining component was 
recorded, it was not recoverable from the magnetic tape. 
Hicmecasncmml or Lats Malfunction was che failure of the 


mer OsvteGiec = wiiten could not be replaced during the 
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experiment. Availability of the vertical velocity data 

would have enabled a cross-correlation investigation of 

PmcomcirTecus Of Particle velocity on the acoustic fluctuations, 

although time for this analysis may not have been available. 
tCmewiine Wile COltOsu lt TOlmcunle  nola@crack current mebers 

are being analyzed in New Zealand and the data corresponding 

to the runs analyzed herein were not available during this 

effort. Hence no information concerning advective processes 

affecting the range could be determined. 

Surface wave data were restricted to SueTiitican tewavic 
height as recorded on the logs. Data on wave periods present 
during a run would have aided the analysis of periodicities 
Poaenme= auLocorrelation functions. 

As detailed in Section I.B., the purpose of the experi-~ 
ment was to investigate both the acoustic fluctuations and 
pae Nrcrteosecames feo ures ol) Che medium. 20 “ssi felit Chat the 
parameters listed above would have been useful in an analysis 
of the temperature microstructure aS well as of the acoustic 
propagation. With regard to-temperature analysis, a measure 
of insolation more quantitative than logged comments about 


We wooOmcnurono would have been advanvagceous. 


See kU LENGTH 

ime phase Sie rune Meneth wast nominally stensminuteis. By 
Heine ea maximum time-lagz of one-tenth the run length, periods 
Liriano lO about 55 seconds were able to be detec- 


LeGmwith weasonablew aefinition in the autocorrelation and 





spectral energy density functions. Better resolution could 
have been obtained with the same number of lags by decreasing 
the maximum lag time or increasing the run length. Decreas- 
ing the maximum lag time has the effect of shortening the 
longest resolvable period with the same number of lags. 
Thus, maximum lag was selected as"large as possible to take 
advantage of the entire run length as a fixed parameter. 

During Phase Il, run length was increased to twenty 
minutes to enable the detection of longer periods of fluc- 
tuation - up to two minutes if a maximum lag of one tenth 
the run length criterion were used. 

The long period oscillation apparent in Runs 116 and 
117 was not detected in the analysis because of the short 
length of the run compared to the period of oscillation. 
MNowever, the wse of run lengths long enough to ensure analy-— 
sis of such long periods would have been prohibitive in 
ie Cicme@iCoboims LOraceereguilrencn ts mamas Wolter. conpurer 
analysis requirements. Since no information was available 
beforehand to suggest the presence of long periods, run 
lengths of 20 minutes were best suited to the purposes of 


the experiment. 


C. RECORDING OF DATA 

Two Lape recorders were used to record separate parameters 
Siesimultaneous Gime series,#al system nou beneficial to 
fanalysisc duewto the problem of tape synchronization during 


mMecenSLeucCpV@newe 1120 a fourteen track tape recorder been 
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available, all of the data could have been recorded simul- 
taneously on one machine. As it was, the surface wave data 
was recorded on a different tape from the one used to record 
the pulse height and temperature data. 

im preparations fior digitization of the data, it was 
discovered that on at least one of the dubbed tapes, the 
pulse height data had exceeded the +40% deviation from the 
carrier center frequency for the FM recording band and tape 
speed used during the dubbing process. Attempts to reproduce 
the data on meeane recorder other than the Sangamo were 
Masuccess#ul™ due to the’ bandwidth limitation. Careful 
Costin Of tChewsangamo tape recorder, which was used to dub 
the tapes, revealed that this machine was capable of handling 
signals that exceeded the standard +40% deviation; in fact, 
the discriminators proved to be linear out to at least +602 
deviation from the center frequency. Had this tape recorder 
not been designed to such broad specifications, much, and 
perhaps a@ligwor the pulse” height SG would have been 
irretrievable due to an apparent Bperacor overs lene i 


recording the data out of band. 


D. PHASE II RUNS 

ics temberavure Microstructure activity during the runs 
conducted in Phase II was generally higher than during the 
maasc runs. However, ther acoustic Pulse repetition rate 
employed in Phase II precluded an analysis of the pulse 


helene time serles=oy the digital methods used during this 





thesis investigation. With a repetition rate of 3 Hz, a 
pulse was transmitted every 333 milliseconds. Nominal travel 
time from the projector to the hydrophones was 195 msec, after 
which a pulse height was held by the sample and hold circuit 
for, 200 imscce Lh, theetrigger signal that initiated ga 
pulse sequence arrived at the sample and hold circuit while 
the previous pulse height was still being held. The result 
was a distorted, two-step, pulse shape with the actual pulse 
height not readily distinguishable from the effect of The 
trigger arrival. A minimum pulse repetition interval not 
less than the sum of the hold time (200 msec) and maximum 


travel time would have precluded this problen. 


E. REFLECTED SIGNALS 

Calculated reflected path arrival times varied from 
0.2 to 0.8 msec after the direct path arrival for the various 
PCmuLC CUCM paLNSseem FOr the smalleSpatime ditference, O.c 
Misee, ony { cyclesiot the 35kHz direct path pulsemwould 
have arrived before the reflected signal reached the hy dro- 
phone. Due to these considerations, the sample and hold 
circulby Was designed to require only fivescycles) to sense 
SCO lCmemom Si CicimmnuS~ Wilton the scirculpecparating 
Deol: Ween Lerierencemerom ret lectedgsignals# should 
Pagomoccwm@oCrmmmnOvwe Vom surtace refilecred signals were 
observed on Hydrophone 1 under certain circumstances (low 
tide and 35 kHz frequency). This problem was alleviated 


Bye reuuce@on Jn mast helen. 
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F, LOG KEEPING 

Reconstruction of an experiment that was conducted one 
year prior to analysis depended heavily on the comments 
meecorded on the logs.  Scrupulous attention to detail in 
original log-keeping made this effort relatively successful. 
However, certain omissions occurred. These included lack 
of weather data for a run, no notation of the acoustic 
frequency used during a run, and inadvertent omission of 
Once Ful Where uransterrangs dava from the rough log to smooth 
lee. sihis qa the human factor in experimental work and 
only constant care and cross-checking can minimize errors 


such as these. 





V. SUMMARY, RECOMMENDATIONS, CONCLUSIONS 


The experiment conducted was unique for several reasons: 

1) It involved simultaneous transmission along near-ortho- 
Sonal pathssevto investigate the degree of isotropy in acoustic 
propagation. 

2) Two acoustic frequencies were employed to investigate 
the frequency dependence of the acoustic fluctuations. 

3) A unique sample-and-hold device was designed and con- 
structed to measure acoustic pulse heights without the need 
for envelope detectors or recording of the pulse train. 

4) A unique micro-scale temperature sensor system was 
also designed and constructed to measure thermal microstruc- 
bag 

5) The sensors were fixed to stable platforms, and were 
remotely operated and monitored. 

6) Multiple sensors were employed to investigate all 
important aspects of the ocean environment. 

7) Because of the fixed locations, experimental runs 
were able to be conducted at various times of the year at 
PACA Same 5 LL.e.. 

iieCeapDaGtuEGiemvncCme YDCM Men ana tlwied im bas thesis repre— 
ScibeadnOnMle mene anivial Sivages Of a project that has since 
expended sddagsa1s Still in operation. The fact that the pro- 
meCt Was Sulu) 20 16S embryonic development when these runs 


were conducted resulted in several shortcomings which have 
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already been commented upon. However, it is felt that useful 
information about the environment and its relation to the 
acoustic propagation results presented herein can be gained 
from implementation of the following recommendations: 

1. Conduct an analysis of the temperature microstructure 
data for the same runs as analyzed here. Comparison of the 
results of spectral analysis may indicate the degree of 
influence of surface wave activity on the microstructure. 

2. Obtain the magnetic tape records for the surface wave 
spectrum for ecmeley: analysis and comparison to the results 
of acoustic and temperature analysis. Because of previously 
moted problemsijwith tape synchronization, this effort would 
Beeoaply not prove conclusive. 

For future experiments with purposes similar to this 
ewe Ue los Suesesited tira ; 

1. Every effort be made to measure turbulent velocity, 
temperature microstructure, surface waves, and acoustic 
fluctuations simultaneously. It fe imperative that all 
sensors be maintained to accurately define the environmental 
iia Le are TevOIms « 

2. All data be recorded on the same magnetic tape simul- 
taneously to avoid problems with synchronization. 

3. Vertical and horizontal thermistor arrays be employed 
vice a single thermistor. 

i. The geometry of the acoustic range be carefully 
Mivestaecaypauewiuhn trelard to the efTect of reflected path 


Signals. 
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5. Equipment operating characteristics be well understood 
to avoid difficulties such as those experienced in this 
experiment with tape recording technique. 

6. Vertical temperature profiles be included as a data 
requirement. 

7. Run length be modified as necessary to fully observe 
the longest period fluctuations suspected to be present. 

Conclusions from this analysis of the acoustic variabil- 
ity are listed below: 

lL. The awenik fluctuations were caused by variations 
in the refractive index along with each propagation path. 
These variations were most likely caused by the temperature 
mreroscructure. and by turbullent water particle velocity 
resulting from surface waves. Whether the turbulent velecity 
served solely as a mechanism for temperature distribution 
aodamoditacapton, or whether if also caused changes in 
moiee rel racvivesandex by advecting other inhomogeneities) 
peen as) bubbles or sediment part lene COMMNG. vais |e cline 
Hinedunron thissanalysis. 

e. The fluctuations exhibited various degrees of anisotropy 
With the propagation along the acoustic path oriented perpen- 
dicular to the predominant swell direction generally subject 
momesveateravariability. 

Se. No reliable trend wais obServed to indicate a Significant 
mete rence an Variability due solely to the acoustic fre- 


quencies used. 
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4. The possible influence of internal wave activity was 
indicated by a long period fluctuation of the mean in at 


least two runs. 





APPENDIX A 


Pulse Height Time Series and Analysis Results 


This appendix contains two sections of figures. The 
first section presents representative samples of the tempera- 
ture microstructure activity present during each run. On 
each sample, time increases to the right. Appropriate 
scales and the qualitative classification are noted for 
each figure. 

The second section presents four plots for each of the 
twelve runs analyzed. With two exceptions, the X and Y 
scales for the same plot from different runs are identical 
to ald comparison among runs. The first exception is for 
mhUMN sj sewiteh WaS twice as long as any other; for this run 
the X-scale of the time series and autocorrelation plots 
has been modified and is defined on the plots. The other 
exception applies to the Y-scale of the spectral energy 
density plots. This scale varies slightly among the pilots 
due to the varying degree of fluctuations. The density 
function values are acceptable for comparison among runs. 

Reference to Tables I - IV is helpful in obtaining a 


meaningful comparison among runs. 
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